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ABSTRACT 

We investigate the stability properties of a hot, dilute and differentially rotating weakly 
magnetized plasma which is believed to be found in the interstellar medium of galaxies 
and protogalaxies and in the low-density accretion flows around some giant black holes 
like the one in the Galactic center. In the linear MHD regime, we consider the combined 
effects of gyroviscosity and parallel viscosity on the magnetorotational instability. The 
helical magnetic field is considered in the investigation. We show that the gyroviscous 
effect and the pitch angles cause a powerful gyroviscous instability. Furthermore, in 
most of the cases, plasma with the above mentioned properties is unstable and the 
growth rates of the unstable modes are larger than that of the magnetorotational 
instability. 
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1 INTRODUCTION 

Enhanced outward angular momentum transport problem 
in accretion discs have been treated in hydrodynamic (HD) 
and magnetohydrodynamic (MHD) contexts by referring to 
the fluid description as well as the kinetic approximation. 
An overview of the past efforts on solving the problem is 
given by Balbus & Hawley (1998) and Balbus (2003). 

Shakura and Sunyaev (1973) appealed to the large 
Reynolds number to get over the difficulties presented by in- 
adequate kinematic viscosity. The breakthrough eventually 
was made by Balbus and Hawley (1991) when they showed 
that differentially rotating an accretion disc with angular 
velocity decreasing outward and threaded by a weak mag- 
netic field is linearly unstable. This instability is known as 
the magnetorotational instability (MRI). Balbus and Haw- 
ley (1998) showed that even non-Keplerian "thick" discs are 
unstable in the presence of a weak magnetic field and thus 
claimed the more generality of the MRI. The literature con- 
cerning the applicability of the MRI in linear as well as 
nonlinear regimes on various astrophysical systems is rich 
in content. The reader may refer to the review article by 
Balbus (2003). 

Balbus (2004) investigated the instability of the differ- 
entially rotating dilute plasma and showed that the maxi- 
mum growth rate of the magnetoviscous instability is greater 
than that of the MRI. Also, he noted that the parallel vis- 
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cous stress can lead to accretion disk turbulence, even if 
Lorentz force is negligible. Balbus refers to the resulting vis- 
cous tensor, in the limit of LUdTi ^ 1 , as the Braginskii 
viscosity which he claims to be appropriate for interstellar, 
galactic and protogalactic discs. But Balbus (2004) did not 
take gyroviscous force which is one of the components of the 
stress tensor and due to finite Larmor radius (FLR) effect 
into account. 

Quataert, Borland and Hammett (2002) explored the 
kinetic version of the MRI which is the regime where the 
wavelengths are much larger than the proton Larmor radius 
in a hot accretion flows onto compact objects. They showed 
that in the various limits of plasma beta, growth rates dif- 
fer considerably. The major conclusion they draw is that 
kinetic effects they considered do not change the stability 
criterion from the MHD resuh, but do change the growth 
rate significantly at high beta. 

Ramos (2003) presented an analysis on the dynamic 
evolution of the parallel heat fluxes in a collisionless mag- 
netized plasma. Fluid description of collisionless plasmas 
was first given by Chew, Goldberger and Low (1956). They 
worked in the lowest order or zero Larmor radius limit. Dia- 
magnetism and other multifluid effects require higher or- 
ders in the gyroradius expansion. Taking into account the 
above mentioned shortcomings facing the collisionless mag- 
netised plasmas, Ramos (2005) derived the first-significant- 
order FLR systems of fiuid moment equations. He claims 
that this formahsm can account for the gyroviscous stress, 
the pressure anisotropy and the anisotropic heat fluxes in a 
plasma with an arbitrary magnetic field geometry. 
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Islam and Balbus (2005) generalized the viscous insta- 
bility examined in Balbus's earlier study (Balbus 2004) by 
including the dynamical effects of magnetic tension force. 
They found that the growth rates are lower than the case 
when tension force is not considered. Nevertheless, they ar- 
gue, the growth rates are still higher than the maximum of 
the standart MRI. They conclude that the magnetoviscous 
instability may find an application in galactic discs and halos 
and low density accretion flows around the Galactic center. 

Fcrraro (2007) considered the hitherto neglected ion gy- 
roviscosity effect which represents the first-order FLR cor- 
rections to the two-fluid MHD equations. He showed that 
FLR effects are much more important than the Hall effect in 
the limit of weak magnetic fields. Also, Ferraro claimed that 
in the collisional limit, gyroviscous effects at scales much 
larger than tli may completely stabilizes the MRI. But he 
considered the magnetic field configuration with Bo = Bz 
that the parallel viscosity doesn't play a role in the MRI. 
Although it is clear that extending his analysis to a more 
general magnetic field configuration (especially that the par- 
allel viscosity is important) complicates the analysis consid- 
erably, it is necessary to understand the true nature of this 
instability. Also, FLR effects that produces coUisionless vis- 
cosity represent "non-ideal" MHD terms like Hall term in 
the fluid formalism. Therefore, the FLR effects may intro- 
duce important corrections to the dispersion relations for 
MHD instabilities and waves. An investigation considering 
both parallel viscosity and gyroviscous effects is the primary 
motivation of this paper. 

The plan of the paper is as follows: in Section 2, we 
give our two fluid equations containing hithereto unexplored 
hot, dilute differentially rotating plasma medium threaded 
by a weak (subthcrmal, i.e. magnetic energy density is small 
compared with the thermal energy density) helical magnetic 
fleld wherein the gyroviscosity and parallel viscosity may 
play important roles in stability. In Section 3 we present the 
linearized MHD equations and their solutions and finally in 
Section 4 discussion and conclusion are given. 



2 TWO FLUID EQUATIONS 

In a dilute plasma, ion gyroradius is much smaller than 
the ion collision mean free path, i.e., tli <C Ai, and also 
ion cyclotron frequency (cjci) greatly exceeds ion- ion colli- 
sion frequency {vi), i.e., e = LOciTi ^ 1 where = 1/vi. 
Conventional MHD theory assumes that plasma pressure is 
isotropic for simplicity. However, the one that want to make 
a more accurate model of plasma must take into consid- 
eration anisotropy in the pressure. The anisotropy in the 
pressure can be due to the presence of a strong magnetic 
field under laboratory conditions, but under the astrophys- 
ical conditions, the plasma may have a tcnsorial character 
even in the presence of a weak magnetic field, the result- 
ing anisotropy being determined by the Larmor frequencies 
and the macroscopic velocity gradients of the constituents 
of the plasma. FLR effects are accounted for by the gyrovis- 
cosity, which arc terms in the pressure tensor independent 
of collision frequency (Roberts & Taylor 1962). 

Spitzer (1962) gives e parameter for a hydrogenic 
plasma as below: 



1.09 X 10" ^ Tf'Bu,G 
In A 



(1) 



where n is the proton density in cm~ , T4, the temperature 
in units of 10* K, B^a is the magnetic field in microgauss 
and In A is the Coulomb logarithm. For n < 1 and Ti>^l, 
the condition of e S> 1 is fulfilled oven in the presence of 
a very weak field. Also Islam & Balbus (2005) mentioned 
the one which assumes large Reynolds number will obtain 
that the collision frequency is much larger than an orbital 
frequency but much smaller than the cyclotron frequency 
(i.e., f7 < i^i < 

Under this condition, MHD equations describing the 
plasma dynamics should include the anisotropic transport 
terms arising from the free fiow of the particles along the 
magnetic field lines (Braginskii 1965). The parallel vis- 
cosity by ions is higher than the electrons by the factor 
{mi/meY^^ ■ In order to take into account the arbitrarily 
high order FLR effects the MHD equations should contain 
coUisionless viscous tensor. Balbus (2004) stressed on the 
importance of ion viscosity in a rotating systems. In an at- 
temp to investigate the FLR effect and parallel viscosity in a 
dilute plasma one should consider two fluid equations. The 
equilibrium state is a differentially rotating dilute plasma. 

The two-fluid moment equations together with Faraday 
and Ampere laws, respectively are given as below (Braginskii 
1965): 



It 



ris-^ = -VPs-V-Tls+qaTis (^E -|- 
dt \ 



V X B 



9B 

dt 



-cV X E 



(2) 

j -msUsgiS) 

(4) 
(5) 



J = —V X B = ene(vi - Ve) 

where the subscript "s" stands for electrons and ions. 
The symbols have their usual meanings, d/dt — d/dt+{v-\7) 
is the Lagrangian derivative. 

If we assume quasi- neutrality, i.e., rie = Zm then the 
mass continuity equation is written only for ions. Let us put 
P = Pi + Pe,TI = Tli + lie and assume that me/rm ~ 0. 
With these substitutions the MHD equations become. 



dp 



+ pV ■ 



: 



dB 
'dt 



V X 



„ J X B cVPe cV • He 

V X B 1 1 



(6) 
(7) 
(8) 



Stress tensor 11 = n" -|- 11-^ -|- 11"" is described as the 
sum total of the parallel (||), perpendicular (_L) and the gy- 
roviscous (gv) components (Braginskii, 1965). In a dilute 
plasma perpendicular viscosity is smaller than the parallel 
viscosity by a factor (r^/A)^ where is the Larmor ra- 
dius and A is the mean free path of the particles. 11"" is a 
dissipationless stress which represents the lowest-order FLR 
correction to the fluid ccjuations and measures the changes of 
particle drift velocities across a gyroorbit (Ferraro 2007). We 
use the Braginskii closure for the stress tensor 11 = 11" -l-II"" 
in our investigation (see, e.g., Ferraro 2007). 



Finite Larmor Radius Effects on Dilute Plasmas 3 



0.96 



2vi 

Pi 



(l - 3bb) (b • W • b) (9) 



b X W • (I + 3bb) + [b X W • (I + 3bb)] ' 



(10) 



Here b = B/B, iOd = eB/niiC are the unit vector along 
the magnetic field and the cyclotron firequency, respectively; 
Vi ion collision frequency. W = Vv + (Vv)^ — 2/31 (V • v) 
is the rate of strain tensor. 



K 1 



+Vpar2D 



dn 1 



dlnRu) 

a In it w 
dfi ^2 1 



cosOSvR + cos05v^ + sin 95vz 

1.2 "1 

Svr 



(13) 



3 LINEARIZED EQUATIONS 

We assume that the Hall effect which is represented by the 
second term on the right hand side of the equation (8) is 
negligible for /3 ^ 1 (Ferraro 2007). This effect is the most 
important one in the low-/3 (/3 is the ratio of the gas pressure 
to magnetic pressure.) plasmas. Therefore will not be con- 
sidered in the present investigation. The third term on the 
right hand side of the equation (8) is the thermodiffusion 
term. Due to their higher masses, ions carry the most of the 
momentum. Therefore the last term on the right hand side 
of the equation (8) is negligible. We work in the Boussinesq 
limit and set V ■ v = 0. The velocity v in the equations 
(6)-(8) is the ion velocity. We seek the solution of the dis- 
persion relation to be derived from the set of equations (6)- 
(8). We work in a cylindrical coordinate system, {R,<j),z). 
Space-time dependence of the axisymmetric perturbations 
is assumed to be of the form, exp(ifc_R + ik^ + ujt). The weak 
magnetic field has a helical shape with = Bo cos 9 and 
Bz = BosinO, where 9 = ta,n~^{Bz/B^) is the angle be- 
tween the magnetic field vector and the (j) axis of the coordi- 
nate system. The pitch angle we define as the one between 
the (f> direction of the cyclindrical coordinate system and the 
magnetic field B, not as the angle between the instantaneous 
velocity vector and the magnetic field as defined in the par- 
ticle orbit theory. Bq is the magnitude of the seed field. The 
reason why we assumed Br = is that the presence of a 
finite Br generates a time dependent (Balbus & Hawley 
1991) which in turn makes the analysis complicated. In the 
equilibrium state of the differentially rotating plasma with 
a Keplerian velocity profile, = RQ{R), plasma pressure 
is assumed to be isotropic. 

The linearized equation of mass continuity is as below: 



kRSvR + kzSvz = 



(11) 



The radial, azimuthal and the axial components of the 

linearized momentum conservation equation are given by the 
equations (12), (13) and (14), respectively, 

ujSvr — 20,dvtb + ikr — 
P 

+-^ikR {B^SB^ + BzSBz) - -^ikzBzSBR 



£rRZ^li^^^'2^SvR + ^sm29Sv^ 
-1-2 |a sin^ 9Svz 



2sine^i(§-l)^^;« 
+{2cos9!^ +E)Svz 
-(2^sin6l-l- A)<5t;^ 



(12) 



+V„. 



5P 1 
LoSvz + ikz h -, — ikzBs&B, 

p A'Kp 

F-^^-SvR + F5v4, + CSv. 



(2'^co,9 + E + 4D{^ 



i2 1 



Svr 



+ (4|a sine + 475^1) Ji;, 



' fci ""^ " ' d In J 

+irtRhG5vz 



■ (14) 



and similarly, the radial, azimuthal and axial compo- 
nents of the linearized magnetic induction equation are given 
by the equations (15), (16) and (17) respectively. 



(jjSBr — ikzBzSvR = 



u6B^ — ikzBzSv^, — 



dQ 
l\nR 



5Br = 



uiSBz — ikzBzSvz = 



(15) 
(16) 
(17) 



where Vpar = 0.96kzPi/2i'ip is the inverse time scale of 

the dissipation due to parallel viscosity; Vgyro ~ klPi/Au)ciP 
is the inverse time scale of the gyroviscous stress, is 
epicyclic frequency. Other constants which depends on pitch 
angle 6* are A = sin 6/(1 - 3 cos 29), B = sm9{l -|- 9cos^ 9 - 
3sin2 6>), C = 2 3111^ 9 {-1 + 3 3111^9), D = 33111^90039, 



2cos6'(l + Ssin-" 



2sin6'D - sin26>, G = 



sin 61(1-1-3 cos 29) and = 3 sin 9 cos^ 9. 

The dispersion relation (DR) derived from the set of 
equations (11)- (17) is given below: 



w'* -|- 030;^ + 02W^ + oiw + ao = 



where 



as = GVparS -j^ 



«2 = 2kWA + f + ^V^yro^\2sy' - Af 

+ Vgyro^'{G + H)+ 4'^Vgyro{2Sy^ " A) 



2 2 2 k\ ~ 2 / 

0/1 — 6Vpar6' k^Vj\ ^2 -|- \^parVgyro3s -j^ 

+ '^VpaM2cD + ^Vgyro&y^DQ' 

1.2 

+ {V9yrof3y^D{2sy' - A)^' ^ 



(18) 



(19) 



(20) 



(21) 



00 = ■^v'V^yro2D''{Q.'f + [klv\ + %yro^{Gl2 i?)] 



klv\ + Ygyro^Gl2 + %^i 



(22) 
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where y = kn/kz, O,' = dO? /dlnR, s — sin9, c — cos9, 
kj_ = + k1 cos^ 6 and Vgyro = klPi/AQuJciP- If we set 11 
the stress tensor to zero and assume = 0, our dispersion 
relation is reduced to the one given by Balbus & Hawley 
(1991) (taking into account the difference in their definition 
of omega; see their equation 2.9). If we set the gyroviscous 
force to zero, then wc recover the dispersion relation given 
by Islam & Balbus (2005) (see their equation 31). 



3.1 Instability Criterion and Numerical Solutions 

Equation (18) describes four low-frequency modes that exist 

in a plasma properties of which described above. A necessary 
and sufficient criterion for instability is clearly given by the 
Routh-Hurwitz theorem 



«0 = J^y^Vsyro^D^n'f + [klv\ + %yron'{G/2 + H)] 



2 

z r^i 



ktvl + Vgyro^'G/2 +YZ^ 



< (23) 



In the absence of gyroviscosity effect, this criterion is 
recovered as (k • va)^ < — fi' and it is an ideal MRI cri- 
terion which states "the combination of a negative angular 
velocity radial gradient with almost any small field will lead 
to dynamical instability" (Balbus & Hawley 1991). When 
we consider FLR effect the instability criterion (23) is very 
complicated. Because it depends on magnetic tension force, 
gyroviscosity force, angular velocity gradient and 0-pitch an- 
gle. Also, angular velocity gradient, gyroviscosity force and 
constants depend on pitch angle may bo positive or negative. 
For the simplicity wo consider only B = Bqz case. In this 
situation pitch angle is 90° and D = 0, G = —2, iif = 0. If we 
assume Q! < 0, the instability depends on the sign of Vgyro- 
For the case of Vgyro > 0, since the first factor of the second 
term is positive, the instability criterion is reduced to the 
condition that the second factor must be negative. Although 
gyroviscous force is coupled with differential rotation which 
is the source of free energy of MRI, it acts in the same direc- 
tion as the magnetic tension force. Therefore it suppresses 
the instability. For the case of Vgyro < 0, the gyroviscous 
force acts in the opposite direction of the magnetic tension 
force and enhances the instability. On the other hand, for 
Q' > 0, the situation is reversed. As a result, the gyrovis- 
cosity force stabilizes or destabilizes according to whether 
(f2 • B) is positive or negative. In this regard, the FLR ef- 
fect is similar to the Hall effect. The contribution of the first 
term (in the inequality (23)) to the instability arises only 
when fei? is taken into account and 0° < 6^ < 90°. This term 
is always positive, therefore it has a stabilizing effect. When 
wc consider different pitch angles and wavenumbers, the in- 
stability criterion is highly complicated due to the change of 
the signs and relative magnitudes of the terms. Therefore, in 
order to reach a definitive results the numerical solutions of 
the dispersion relation (18) are required. We have rewrittten 
the dispersion relation in its dimensionless form, i.e. all the 
terms of Eqn.l8 is divided by O,*: 

(24) 



7" + &37^ + &27^ -I- 617 -I- 60 = 

where 



(25) 



.dlnn^ 
'avro^ dlnJi 



&2 = KyroX^^^^^iG + H) + 4^Vg''yroX{2sy^ - A) 



+2X + ■^k' + ^{Vglrofx\2sy' - Af (26) 



61 = 6V"s^ktX -i + V"VZr„X''3s''G 



dlnR A;2 

k'i - ^ dlnfl\ ^ , ki~^„ a^Yn^h^ 



y par'' 
1,2 J i„ r>2 



par * gyro-' 
,2 



+{V^yrofX^iyD{2sy' - A) 



d\nR P 



(27) 



60 = 



^ 1 + Ky 



dlnR 



G/2 + 



fc2 dlnn^ 



fc2 dlnR 



Xll + Kyro^{G/2 + H) 



kl 2,,~.n x2^2^r^2/rflnn\2 



V^y\KvrofX^2D\ 



-r 



(28) 



yn 
*' par 



jfc2" v™-"/ V (^inii 

where 7 = ui/Q., X = kzv\/Q^, Vgyro = Vgyro /X and 

Vpar / ^X . 

Vpar ^ ^g"r-o relationship is easily derived from con- 
dition of dilute plasma e = oJdTi 3> 1 where Ti = 
Taking into consideration this relationship, all the figures 
are drawn as dimensionless growth rate versus wavenumber 
for Vpar = 1000 and for the Keplerian rotational profile. The 
growth rate depends sensitively on the pitch angle and the 
gyroviscosity parameter Vgyro/^X = V^yro- 

Figure la shows the fastest growing mode (i.e., kn = 0) 
in the case when the rotation axis (fJ) and the B^ are par- 
allel (f2 tt Bz); and Figure lb is plotted for kn/kz = 1 and 
also f2 tt Bz. Figure la and lb are drawn for the V^yro = 1- 
Although all the wavenumbers arc unstable for the pitch an- 
gles 15° and 45° when kn = and kn/kz = 1, for the 9 > 60° 
all the wavenumbers are stable only when kn/kz = 1. Also, 
the growth rates exceed the maximum growth rate of MRI 
(0.75n) especially for the 6 = 30° in the both cases. When 
the pitch angle between the magnetic field lines and the 
angular velocity gradient is 90° there is the minimum vis- 
cous transport of angular momentum (Fig. la). Instability 
tends to grow when the angular velocity gradient acquires 
a component along the magnetic field lines. But for greater 
pitch angles (> 45°) the radius of curvature becomes greater, 
therefore tension force becomes smaller. In a dilute plasma, 
the dynamic effects caused by the magnetic tension force 
would lower the growth rates (Islam & Balbus 2005). One 
would expect greater growth rates when the tension force 
is smaller. But, as Islam & Balbus(2005) states that the 
dependence of the growth rate on the tension force (in our 
opinion, possibly through the pitch angle) shows the compli- 
cated nature of the magnetic stresses "as both a stabilizing 
and a destabilizing agent" . Therefore it seems as though the 
viscous angular momentum is more effectively transported 
when the pitch angle is within the range of about 15° — 45° . 

When the angular velocity vector and the Bz compo- 
nent of the magnetic field are anti parallel (fl t-i B^), gy- 
roviscosity assumes negative values. Figure Ic and Id are 
drawn only for the V^yro = — 1 in the cases of fcii = 
and ka/kz = 1, respectively. In the case of = all 
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Figure 1. The growth rate of gyroviscous instability for the cases of f2 '[^ Bz (Fig. la and lb) and f2 '['J, (Fig. Ic and Id). Figure 
la and Ic is drawn for kji = and Figure lb and Id for kji/k^ = 1. In the situation V^y^^ = 1 for the pitch angles are smaller than 60° 
and in the situation V^^^ = —1 for the pitch angles are greater than 60° the instability spreads over a wider range of wavenumbers. 



the wavenumbers are unstable only for the 8 > 45° (Fig. 
Ic), in the case of kn/k^ = 1, they are unstable for the 
9 > 60° (Fig. Id). In addition, in Figure Ic and Id maxi- 
umum growth rates are larger than the ideal MRI growth 
rates for the great values of pitch angle. Besides, when we 
consider finite kn wavevector, an unstable mode emerges for 
e = 90° 

The normalized wavenumber of the fastest growing 
mode is given as k^VA/^ = (15/16)^''^ (Balbus & Hawley 
1992). We adopt the same value in the present investigation. 
We plotted dimensionless growth rates versus pitch angle 
and kn/kz in Figure 2. For the V^y^o = 1, growth rates are 
given in Figure 2a and for the V^yro ~ they are given in 
Figure 2b. In the case of ft growth rates have max- 

imum values for the small pitch angles and kn values (Fig 
2a) . While there is no unstable wavenumber for the 9 > 60° 
(Fig 2a) for the case of f2 tJ- these pitch angles have 
maximum growth rates (Fig 2b). 

To see the effect of the gyroviscosity parameter on the 
growth rates we assumed the pitch angle as 9 — 45°. Then, 
we plotted dimensionless growth rates versus gyroviscosity 
parameter and wavenumbers for the kn/kz = 0, kn/k^ — 1 
and kzVA/ii = (15/16)^''^ respectively in the Figure 3a, 3b 
and 3c. We clearly see that all wavenumbers are unstable for 
the positive values of the groviscosity parameter, i.e. in the 



case of n tt Bz. The positive values of the gyroviscosity 
parameter have maximum growth rates case in the large 
wavenumbers, but again the growth rates are greater than 
0.75!^. 

In Figure 4, we assumed the pitch angle as ^ = 60°. 
Again, we plotted dimensionless growth rates versus gyro- 
viscosity parameter and wavenumbers for the kn/kz = 0, 
kn/kz — 1 and kzVA/^ = (15/16)^''^ respectively in the 
Figure 4a, 4b and 4c. We clearly see that all wavenumbers 
are unstable for the negative values of the groviscosity pa- 
rameter, i.e. in the case of f2 Bz. Maximum growth rates 
emerge in the small wavenumbers and and the negative val- 
ues of the gyroviscosity parameter contrary to previous case 
(Figure 3). 



4 CONCLUSIONS 

In this paper we have investigated a linear axisymmetric 
analysis of the MRI with FLR effect and paralel viscosity in 
the hot, dilute and differentially rotating weakly magnetized 
plasma. In a dilute plasma, ion cyclotron frequency much 
exceeds the ion-ion collision frequency. In this regime, the 
viscous stress tensor is given as the total of the parallel and 
the gyroviscous components (Braginskii 1965). This work 
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Figure 2. The dimensionlcss growth rates of instability versus pitch angles and kn/k^ for the normahzed wavenumber of the fastest 
growing mode, i.e. k^VA/^ = (15/16)^/^. In the situation V^y^o = 1 (Fig- 2a) for the pitch angles are smaller than 60° and in the 
situation V^y^o — ~^ (Fig. 2b) for the pitch angles are greater than 60° all the wavenumbers are unstable. 




kv./n kv./n kJk 

A A R z 

(a) (b) (c) 

Figure 3. The dimensionless growth rates of instability versus gyroviscosity parameter and wavenumbers for the fc/j/fcz = 0, fc^/fc^ = 1 
and fc^D^/n = (15/16)^'^^ respectively in the Figure 3a, 3b and 3c. The pitch angle is taken 45° in these figures. For the positive values 
of gyroviscosity parameter, i.e. in the case of f2 Bz, all the wavenumbers are unstable. 



is not only extension of the Islam & Balbus's (2005) study 
which is taken into account only parallel viscosity but also is 
the extension of the Ferraro's (2007) study which is included 
parallel viscosity and gyroviscosity only with B = Bz mag- 
netic field geometry (also parallel viscosity is negligible in 
this geometry). We investigated the nature of gyroviscous 
instability with the more general magnetic field configura- 
tion, i.e., a helical field for which the parallel viscosity is 
important. We obtained the dispersion relation of this in- 
stability and also derived the instability criterion. 

The gyroviscous instability shows a complicated depen- 
dence on the geometry of the magnetic field through the 
pitch angle 9. Quataert, Borland and Hammett (2002) also 
draw attention to the dependence of the growth rates on 
the orientation of the magnetic field and the wavevector of 



the mode. More important than that is the complex coupling 
between the dynamical effects and the geometry of the mag- 
netic field that makes it almost impossible to single out the 
net effect of any one agent contributing to stability or insta- 
bility. But one thing is certain that in its more general case 
the magnetized, hot and differentially rotating plasmas are 
always unstable. 

Our results show that when only gyroviscosity act on 
the hot dilute and differentially rotating plasma it brings 
about instability and growth rates of the instability are 
higher than MRI (see Fig 1 and 2). In this situation, the 
finite kn is destabilized all the wavenumber interval for 30° 
in the case of ft tt (Fig. lb and 2a) and for pitch an- 
gles greater than 60° in the case of ft (Fig Id and 
2b). Although parallel viscosity is greater than gyroviscosity, 
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Figure 4. The dimensionless growth rates of instabiUty versus gyroviscosity parameter and wavenumbers for the kn/kz = 0, kn/kz = 1 
and kzVj^/Q = (15/16)^/'^ respectively in the Figure 4a, 4b and 4c. The pitch angle is taken 60° in these figures. For the negative values 
of gyroviscosity parameter, i.e. in the case of SI fj, B^, all the wavenumbers are unstable. 



we see that parallel viscosity slightly supresses the unstable 
modes for the cases of ft Bz and ft ti according to 
the situation which parallel viscosity is neglected. Also there 
is an interesting situation: When the angular velocity vec- 
tor and the Bz component of the magnetic field are parallel 
(n tt B2), pitch angles smaller than 60° have maximum 
growth rates in the great wavenumbers; when the angular 
velocity vector and the Bz component of the magnetic field 
are anti parallel (ft B^), pitch angles greater than 60° 
have maximum growth rates in the small wavenumbers (see 
Fig. 3 and 4). 
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In summary, as shown in the figures, the gyroviscos- 
ity effect and pitch angles cause a powerful instability. This 
gyroviscous instability has a very large gowth rates with 
O.Sn — for the different pitch angles. The gyroviscosity 
instability may be candidate for amplifying very small seed 
fields. Therefore, this instability may apply to protogalax- 
ies, low-density accretion fiows and the intracluster medium 
of galaxy clusters which the magnetic fields are sufficiently 
small so that dilute plasma condition is satisfied. 

The last but not in the least the dependence of the 
growth rates on the pitch angle forces us into a qualititive 
discussion. We assumed helical magnetic field geometry. E x 
B drift generates no current, unless we consider an orbit- 
averaged E X B wherein the Larmor radii of electrons and 
ions are different added to this current and viscous fiows 
generates a net current. Depending on the pitch angles of 
current density and the magnetic field vectors the plasma is 
either pinched or diluted locally by the Lorentz force, J x B. 
If the pitch angle of the current density vector is smaller than 
the magnetic field pitch angle the Lorentz force will dilute 
otherwise pinch the plasma locally. It appears that when the 
relative pitch angles of J and B makes the growth rate of 
the unstable mode higher for all pitch angles range, "the 
currents caused by FLR effects becomes in phase with the 
current of the ideal-MHD MRI eigenmode" (Ferraro, 2007). 
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